Bounds on the dragging rate and on the rotational 
mass- energy in slowly and differentially rotating 

relativist ic stars 

M. J. Pareja 

Dept. de Ffsica Teorica II, Facultad de Ciencias Ffsicas, 
Universidad Complutense de Madrid, E-28040 Madrid, Spain 
g; e-mail: mjpareja@fis.ucm.es 

O 

(N 

P-c Abstract 

CD ■ 

00 | For relativistic stars rotating slowly and differentially with a positive angular 

lO ' velocity, some properties in relation to the positiveness of the rate of rotational 

dragging and of the angular momentum density are derived. Also, a new proof for 
the bounds on the rotational mass-energy is given. 



> 
00 



X 



PACS numbers: 04.40.Dg, 97.10.Kc, 02.30. Jr 
I. INTRODUCTION 



ON 

o 

CO 

In the prescription for calculating a slowly and differentially rotating relativistic stellar 
configuration the field equations are expanded in powers of a fluid angular velocity pa- 
rameter and the perturbations (around a non-rotating configuration) are calculated by 
(3J[)! retaining only first and second order terms. Hartle^ has derived these equations of struc- 
ture in the rigidly rotating case. However, at first order the only effect of the rotation 
is to drag the inertial frames; at second order it also deforms the star. In fact, the first 
order equations of structure reduce to the time-angle field equation component (to first 
order), which is a partial differential equation, linear in the angular velocity of cumulative 
dragging (or dragging rate). 

This linearity in the dragging rate potential persuades us to write that equation 
in appropriate coordinates — in order to avoid the coordinate singularity occurring on 
the axis of rotation in spherical polar coordinates, generally used in the slow rotation 
approximation — so that in the new coordinates the equation writes in a "regular" form 
as an elliptic equation with bounded coefficients, and to apply a minimum principle for 
generalized supersolutions in the whole domain (interior and exterior of the star). Mak- 
ing use of the asymptotic flatness condition, this will lead us directly to the positivity of 
the dragging rate, provided that the distribution of angular velocity of the fluid is non- 
negative everywhere (and non-trivial) and that we start from a reasonable unperturbed 
(non-rotating) stellar model satisfying the weak energy condition. This and the positive- 



1 



ness of other quantities also linear in the angular velocity, like the angular momentum 
density, will be the purpose of this paper. 

The rotational mass-energy, derived by HartleP although accurate to second order 
in the angular velocity, involves only quantities which can be calculated from the first 
order structure equation (time-angle component of the Einstein equations) as well. A 
proof of the positivity and an upper bound on this rotational energy is given in the same 
paper p however using an expansion in eigenfunctions. We shall give here a much simpler 
proof of these bounds, without using that expansion, and, hence, avoiding the non-trivial 
mathematical problems on the existence of these eigenfunctions. 

The paper is organized as follows. After a description of the relativistic rotating 
stellar model in Sec. II, and a brief revision of the concepts of angular momentum density 
and rate of rotational dragging in Sec. Ill, in Sec. IV we concentrate on the slow rotation 
approximation, particularly on the first order perturbations of the metric (linear correction 
of the dragging rate, with description of the unperturbed (zero order) configuration), 
and explicit expressions for the expansions of the angular momentum density and of 
the rotational mass-energy are derived. In the same section the null contribution (at 
first order in the angular velocity) of the integrability condition of the Euler equation 
is discussed, and the time-angle component of the Einstein equations (to first order) 
is written in appropriate coordinates, as a background allowing to apply a minimum 
principle and obtain the first of the properties mentioned above and proved in Sec. V, 
and consequences of that one. Apart from this, as an independent result, an alternative 
proof of the bounds on the rotational energy is given. Finally, in Sec. VI, concluding 
remarks are briefly stated. 

II. THE RELATIVISTIC ROTATING STELLAR MODEL 

The spacetime of a rotating relativistic star is represented by a Lorentzian 4-manifold 
(A4,g) which satisfies the following 

A. Assumptions 

i. the spacetime is stationary in time and axially symmetric, which means that g 
admits two global Killing vector fields, a time-like future-directed one, and a 
space-like one, with closed trajectories, rj, except on a time-like 2-surface (defining 
the axis of rotation) where 77 vanishes; 

ii. the spacetime is asymptotically flat; in particular, g(£,£) — > — 1, g^,??) — > +00, 
and g(£, rj) — > at spatial infinity (the signature of the metric g being ( — h ++)); 

iii. the matter — confined in a compact region in the space (interior), with vacuum on 
the outside, so that (ii) holds — is perfect fluid, and therefore the energy-momentum 
tensor (source of the Einstein equations) is written as 

T = (e + p)u b <g>u b + pg, 

where e and p denote the energy density and the pressure of the fluid, respectively; 
and u b denotes the 1-form equivalent to the 4- velocity of the fluid u (in the exterior 
T = 0; hence, e + p = p = there); 
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iv. the fluid velocity is azimuthal (non-convective) (circularity condition) , i.e. 

u b A £ b A rf = ; 

v. (A4, g) satisfies Einstein's field equations G = 87rT for the energy- momentum tensor 
T of a perfect fluid (hi), where G = Ric — 1/2 Rg denotes the Einstein tensor - 
equations which can also be written in the form 

Ric = 87r(T-itr(T)g); (1) 

vi. e and p satisfy a barotropic (one-parameter) equation of state, e = e(p); 

vii. £ + p > (weak energy condition for perfect fluid, assuming e > 0)P 

viii. the metric functions are essentially bounded. 

B. Form of the metric 

Assumptions (i) and (ii) imply that the two Killing fields commute, [£, 77] = 0,™ which 
is equivalent to the existence of coordinates x° = t and x l = <p such that £ = d t and 
77 = 8$; moreover, by the circularity condition (iv), the spacetime g admits 2-surfaces 
orthogonal to the group orbits of the Killing fields (orthogonal transitivity).™ We may 
then choose the two remaining coordinates (x 2 ,x 3 ) in one of these 2-surfaces and carry 
them to the whole spacetime along the integral curves of £ and 77; accordingly, the metric 
can be written in the form 

ds 2 = g a pdx a dx 13 = 

= gudt 2 + 2g t< f,dtd(j) + g^dcf 2 + g22(dx 2 ) 2 + 2g 2 3dx 2 dx 3 + g 33 (dx 3 ) 2 , 

where the metric coefficients are independent of the time x° = t and azimuthal x 1 = (f) 
coordinates corresponding to the two Killing fields; that is, g a p = g a/ 3(x 2 ,x 3 ). 

When solving Einstein's field equations it is convenient to specify the coordinates x 2 
and x 3 in such a way as to simplify the task; a particular choice, usually made when study- 
ing slowly rotating configurations P is the one which makes (723 = and #33 = g^, sin~ 2 x 3 . 
Hence x 2 and x 3 are chosen so that at large spatial distances the asymptotically flat met- 
ric is expressed in terms of spherical polar coordinates in the usual way. In the resulting 
coordinate system, with the notation x 2 = r, x 3 = 9, and with new symbols for the metric 
functions, the line element reads 

ds 2 = -H 2 dt 2 + Q 2 dr 2 + r 2 K 2 [d6 2 + sin 2 # {d<j> - Adt) 2 ] , (2) 

where H, Q, K and A are functions of r and 9 alone. In these coordinates (r > 0, < 
9 < 7i") the spatial infinity is given by r — >• 00, and the axis of rotation (d^ = r\ = 0) is 
described by 9 — > or n (r > 0). 

Notice that the function A appears in the metric (j2J) as the non-vanishing of the (tip) 
metric component of a rotating configuration. A is actually the dragging rate potential 
(cf. Sec. III). 
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C. Differential rotation 



According to assumption (iv) of Subsec. II. A — circularity condition — , the fluid 4- velocity 
u has the form 

u <t> dd> 

u = u l d t + trcL> = uUdt + fi<9A where = — = — 

w at 

is the angular velocity of the fluid measured in units of t, i.e. as seen by an inertial observer 
at infinity whose proper time is the same as the coordinate t (observer d t ), and w* is the 
normalization factor, such that g(u, u) = — 1, i.e. (w*) -2 = — (g a + 2Qg t ^ + VPg^). The 
star's matter rotates then in the azimuthal direction 0. 

We consider a star rotating differentially, with a prescribed distribution of angular 
velocity 

n = n(x 2 ,x 3 ) = n(r,0), 

an essentially bounded function. However, with the assumptions made (Subsec. II. A), 
the rotation profile of the fluid cannot be freely chosen, this shows up in the following. 
We consider the equation of hydrostatic equilibrium, T at3 .p = (integrability conditions 
of the field equations), particularly, its part orthogonal to the fluid 4- velocity u, i.e. the 
Euler equation, 

dp = -(e +p) a, (3) 
where a is the 4-acceleration of the fluid, a" = V u u, specifically, 

a = — dlntt* + rfu^dVt . (4) 

And the integrability condition of Eq. (J3J), taking into account (vi) of Subsec. II. A, e = 
e(p), is da = 0; following, from (jlj), di^u^) A dVL = 0; in other words, the fluid angular 
velocity, Q, is functionally related to the specific angular momentum times u*, 

u l u^ = :f(0) . (5) 

Nevertheless, as will be seen in Subsec. IV.B, in the slow rotation approximation, at 
first order in the angular velocity, Eq. (jSJ) is no restriction on Q(r, 8). 



III. ANGULAR MOMENTUM DENSITY AND DRAGGING RATE 

The total angular momentum of a rotating relativistic star can be definecP from the 
variational principle for general relativity — for an isolated system which is not radiating 
gravitational waves — , but this is showrPto coincide with the geometrical definition — from 
the asymptotic form of the metric at large space-like distances from the rotating fluid — 
(analog to the ADM mass), which for stationary and axisymmetric (asymptotically flat) 
spacetimes is given by the Komar integral for the angular momentumPa surface integral, 
which, when reformulated using the Gauss theorem and the Einstein equations, converts 
into the volume integral over the interior 

j = J T J v a n p dv (6) 
= j Tfadv = j Tfr-gfftfx, (7) 
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where T is the energy-momentum tensor of perfect fluid, rj is the Killing field corre- 
sponding to the axial symmetry, n" is the unit time-like and future pointing normal to 
the hypersurface of constant t, i.e. n = n t dt, with n< > 0, and dv is the proper volume 
element of the surface t = const., i.e. f x dv = Vol, the volume of the body of the star, 
I = interior of the star (t = const.), g = det(g). The invariantly defined integrand of this 
volume integral (0), Tj^rfnp, is what one would naturally define as angular momentum 
density — coinciding with the standard form in special relativity — , and can be calculated 



= n t (e + pjrfuj, [nt = + Vtg H )} 

= n t (e + p)(u t ) 2 (g t<t> + Vtg H ) 

.2 \ 1/2 

. / ill I < I - - (' 

n t (e + p)(u) g^Vt - A), with n t = H 



9<P<f> 

where A is the metric function (cf. ({2)1 ) such that 

g t<t> = -A g H . (9) 

It is remarkable that, since n t > 0, g^ > 0, and we are assuming the energy condi- 
tion e + p > ((vii) in Subsec. II. A), the sign of the angular momentum density (JHJ) is 
determined by the sign of the difference Q — A. 

The metric function A is indeed the angular velocity of a particle which is dragged 
along in the gravitational field of the star, as seen from a non-rotating observer at spatial 
infinity (dt), so that it has zero angular momentum relative to the axis, p^ = 0, 

d<j> _p<t> _ g^pt _ -g H fd<f>\ d<\> _ 



A is called angular velocity of cumulative dragging (shortly called dragging rate)®*" One 
of the purposes of this work is precisely to find appropriate bounds on the uniformly 
non-negative distribution of angular velocity, Q = Q(x 2 ,x 3 ) > 0, of a slowly differentially 
rotating star, so that Q — A > holds; from where the positivity of the angular momentum 
density (JHJ) follows (Property (c) in Sec. V). 

Observe, in the special relativistic limit (g t< f, — > 0, using coordinates (x 2 , x 3 ) which go at 
spatial infinity to the usual flat coordinates, cf. Subsec. II. B), if the fluid rotates uniformly 
with angular velocity Q positive (negative), then the angular momentum density, Eq. (JHJ), 
is uniformly positive (negative). 



IV. SLOWLY DIFFERENTIALLY ROTATING STARS. FIRST ORDER PER- 
TURBATIONS 

By slow rotation we mean that the absolute value of the angular velocity is much smaller 
than the critical value £l C ru = (M/R 3 ) 1 ^ 2 (taking units c = G = 1), where M is the total 
mass of the unperturbed (non-rotating) configuration, and R, its radius; \Q(x 2 , x 3 )\/Q cr u <C 
1. Thus, stars which rotate slowly can be studied by expanding the Einstein field equa- 
tions for a fully relativistic differentially rotating star in powers of the dimensionless ratio 
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"Tj = (10) 

where |f2 maa; | is the maximum value of |f2(x 2 ,x 3 )| (at the interior of the star). 

A. The metric and the energy density and pressure of the fluid 

We assume that the star is slowly rotating, with angular velocity 

n(r,9) = Q{x 2 ,x 3 ) = 0{fi), 

parameter [i given e.g. by (fTU|) . Because the star (stationary in time and axially sym- 
metric) rotates in the (f> direction ((iv) of Subsec. II. A), a time reversal (t — > —t) would 
change the sense of rotation, as well as an inversion in the direction (</> — > —(f)) would do. 
As a result, the metric coefficients H , Q and K (in (j2J)) and the energy density will not 
change sign under one of these inversions, whereas A will do. Therefore, an expansion of 
H, Q and K, as well as of the energy density, e, and, hence, of the pressure, p, in powers 
of the angular velocity parameter fi will contain only even powers, while an expansion of 
A will have only odd ones. 

Accordingly, at first order in the angular velocity, 0(/j), the metric coefficients, and 
fluid energy density and pressure, are 



H = Ho + Oifj 2 ) 
Q = Q + O(/i 2 ) 
K = K + O(fi 2 ) 



e = e + 0(/i 2 ) 
p = p + O{n 2 ) 



but A = w + 0(/i 3 ), (11) 

where H , Q and K are the coefficients of the unperturbed (non-rotating) configuration, 
and ijj denotes the linear (first order) correction in fi of the dragging rate A, so that, from 
Eq. ©, 



gut = -u) {g H ) Q + 0(fi 3 ) . (12) 

We shall keep here only the effects linear in the angular velocity. The only first order 
O(fi) perturbation brought about by the rotation is the dragging of the inertial frames; 
the star is still spherical, because the "potential functions" which deform the shape of the 
star are 0(/i 2 ). 



1. The (unperturbed) non-rotating configuration 

The starting non-rotating equilibrium configuration is described by the spherically sym- 
metric metric in the Schwarzschild form 
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ds 2 = -e u(r Ut 2 + e x ^dr 2 + r 2 (d6 2 + sin 2 6d(j) 2 ) (13) 
= -H 2 dt 2 + Q 2 dr 2 + r 2 K 2 (d6 2 + sm 2 6d(f) 2 ) (A = 0), 

with A(r), or equivalently, the mass m(r) interior to a given radial coordinate r, given by 

x _ 2m(r) = e . A(r) 
r 

and z/(r), together with the pressure Po(r), and the energy density e (r), solutions of the 
system of equations of general relativistic hydrostatics, which for a non-rotating configu- 
ration are: the equation of hydrostatic equilibrium (Tolman-Oppenheimer-Volkoff equa- 
tion) , 

dpo [e (r) +Po(r)][m(r) + Anr 3 p (r)] 

dr [T) ~ r 2 [l - 2m{r)/r] ' 1 J 

the mass equation, 

^(r)=4vrr%(r), (16) 
dr 

and the source equation for u, 

±(r) = 2 - ^(r) (17) 

dr K ' e (r)+p (r) dr K h K ' 

with the initial boundary conditions 



< po(0) = Po c < oo (central pressure), 
m,(0) = , and 

u(0) = v c (constant fixed by the asymptotic condition at infinity), 

this being the prescription for the interior of the star, that is, inside the fluid, r < R, 
R = radius of the surface of the star, determined by Po{R) = 0. Furthermore, we assume 
Po and eo related to each other by a barotropic equation of state, 

£o = ^o(Po) , (18) 

Po l— *■ £ o(Po) a bounded function on any closed interval, and satisfying the weak energy 
condition 

e +Po>0. (19) 

Observe, from Eqs. (|T5j) and (|T§|). p is a decreasing function from the center, r = 0, to 
the star's surface, r = R; in particular, p > and attains its maximum value po c at the 
center. 

In the exterior (vacuum) the geometry is described by the same line element (|13j) . but 
with the metric function v specified and related to A by 

e Hr) = e -A(r) = ]_ _ 2M w>jR; (2Q) 

r 

where M = m(R) is the star's total mass. 
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Notice that this standard form of the non- rotating metric, (j!3|) . is the limit of zero 
rotation of the general rotating metric in spherical polar coordinates (J5]L 

H^H Q = e u/ \ Q^Q = e x/2 , K K = 1, A -> A = ; 

from where the effect of the rotation can be seen as given by the term d(f) — Adt in the 
place of d(p. 

Note, since, as has been seen above in Subsec. IV. A, at first order in the angular 
velocity there is still no effect on the pressure, and on the energy density, conditions (fTSj) 
and (|19|) for the starting non-rotating configuration are (vi) and (vii) of Subsec. II. A (at 
first order). 



B. Euler equation 

It will be important to note that at first order the integrability condition of the Euler 
equation and, hence, Eq. (jSj) are no restriction on Q(r, 6), that shows up in the following. 
Consider the first integral of the Euler equation (JHJ), namely, 



/ P + - ln[( M *)- 2 ] + / ^(O) dQ = const. , (21) 

Jo £ {P> +P 2 (rfi) Jn 

where Eqs. (@J) and (0) have been used, and Q$ is a given constant (changing the value of 
flo simply modifies the value of the constant on the right hand side). The first term in 
Eq. (|2"T]) is a function of the pressure, which is, to this approximation, a function of r, i.e. 
0(1) with respect to the angular velocity; on the other hand, 

(uY 2 = -(g tt + 2Qg t(j> + ft 2 ^) = H 2 - K 2 r 2 sm 2 6 (SI - A) 2 = 0(1 — (SI — A) 2 ) , 

so the second term is 0((SI — A) 2 ) and, hence, 0(/i 2 ); also, since 

t , h2 ^ K 2 r 2 sin 2 6 (SI- A) ., 

u'u, = (u'f 9H (n - a) = IP _ K ^J Hn _' A)2 = 0(tl - -4) . 

T(Sl) = = 0(SI — A), thus, the third term is 0((SI — A) 2 ) as well, and, hence, 0(/j, 2 ). 
Consequently, to 0(/i), the Euler equation reduces to its static (non- rotating) case, and 
indeed we have presumably already used it to get the starting unperturbed solution. 
Therefore, at this order in the angular velocity, Eq. (JHJ) is no restriction on Sl(r, 9). 



C. The angular momentum density 

Using the definition of to, linear correction of the dragging rate, via the expansion of the 
metric coefficient g^, Eq. (fl2|) . and the metric coefficients of the non-rotating configuration 
(fH?|) . we obtain the expansion for the angular momentum density (JHJ) 



T/ V a n p = n t T; = n t {e + P ){u t ) 2 {g t<t) + ng H ) 

= (n t )o (£q +Po)[(-gtt)oV 1 (-u (flto)o + ^ (fiw)o) + 0{j/) 

= e vl2 (e + p )e- u r 2 sm 2 9 (Q - u) + 0(/i 3 ) 

= (e + po) e~ v/2 r 2 sin 2 fl (SI - u) + 0(/i 3 ) . (22) 
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Thus showing also for the first order rotational perturbation that, since we are assuming 
the energy condition Eo +Po > 0, the sign of the angular momentum density (J2*2*j) to 0(/i) 
is determined by the sign of Q — u. 



D. The rotational mass- energy 

In Ref. 2 Hartle has derived the difference in total mass-energy, M rot , between a slowly 
and differentially rotating relativistic star and a non-rotating star with the same number 
of baryons and the same distribution of entropy, namely, 



M mt = ^^ndj + o(fi 4 



where dJ is the angular momentum of a fluid element in the star (to first order in the 
angular velocity), i.e. from (jJJ), 



dJ = TJ{-g) 1 ^ 



x 



o(m) ; 



taking into account (|22jl and piijl . we obtain an explicit expression for the expansion of 
M rot in powers of the angular velocity parameter //, 



M rot = \ ( dr [ d9 2tt(£ + Pa) r A e {x ~ v)/2 sin 3 # Q(Q - u) + 0(/i 4 ) . (23) 
2 Jo Jo 

E. The time-angle component of the Einstein equation 

The (t<f>) field equation component retaining only first order terms in the angular velocity, 
i.e. from Eq. (JTJ), 

takes the form 

d r [r 4 j(r) d r uj] + ^^<9 e [sin 3 #dH - 16n r 4 k(r)[e (r) + p (r)] [u - fi] = 0, (24) 
sin o 

where we have introduced the abbreviations 

j(r) = e-W^+KOl/a and fc ( r ) = e [A(r)-Kr)]/2 _ (25) 

As outlined in Ref. 1, using the 0-order field equations, (fTlj) . ( fTHj) . (fTH|l and (fTTjl . it follows 

Att r [e {r) + Po {r)} k{r) = -j'(r) , (26) 
(where ' = d/dr) which, substituted into Eq. (j2U), yields 



d r \r 4 j(r)d r u;} + ^ d \sm 3 6 d e u} +4r 3 j'{r) u = 4r 3 j"(r) 0(r, 0) . (27) 

sin 6 

We write this differential equation for the dragging rate in the abbreviated form 

Luj = -^ 2 n, (28) 
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with the linear second order partial differential operator L = L — \I/ 2 , where 

L u := -±—d r [r A j{r) d r u] + - , d e [sin 3 9d e u} and (29) 

^ 2 (r) := = i 67T [ £o ( r ) +po ( r )]^) > Vr>0. (30) 

3{r) 

Equation has been used in (jSHJ), and the sign follows from the assumed energy 
condition (|19|h the functions j and A; (J23J) are always positive. \I/ 2 = in the exterior 
(Vr G [R, oof), where vacuum (e Q = p = 0, cf. (iii) in Subsec. II. A) is considered. 

Specifically, we are only interested in solutions uj = u)(r,$) of Eq. (|2*8Jl in [0, oo[x [0, ir], 
which satisfy the boundary conditions 

uj asymptotically flat ( lim co — 0) , (31) 

uj C 1 -regular on the axis of rotation, (32) 

and a matching condition, namely, to be at least a class C 1 function on the surface of the 
star -which is spherical at first order rotational perturbations- 



class C 1 across r = R. (33) 



Notice, (|3*Tj) follows from our star model (condition (ii) in Subsec. II. A), and it can be 
easily seen that ()33)1 follows from the equation itself, provided that 9) and Q(.,8) are 
at least essentially bounded (e L°°) — as has been assumed — i.e. even if they have a 
jump discontinuity. 

At the star's surface, r = R, higher regularity of uj{.,6) is not guaranteed by the 
equation, due to a jump discontinuity of the function \1/ 2 at this point. For this reason, 
we shall be considering (in the following section) generalized (g W 1,2 ) solutions uj of 
Eq. (|28|) in the whole domain (interior and exterior). 

"Coordinate change". 

In order to avoid the coordinate singularity occurring on the axis in polar coordinates 
(r,0), and wishing instead to have in the differential operator (|29jl a Lapacian in some 
higher dimension, we consider the following "change of coordinates" . Firstly, we introduce 
isotropic cylindrical coordinates in the meridian plane, 

(r,0) ^ (p := h(r)sin6 , z := h(r)cos6) GKjxR, (34) 

with the function h satisfying the following ordinary differential equation of first order 
with separated coefficients 

h'{r) _ e^ 2 
h(r) ~ r 

(which makes the coefficient of the crossed derivatives in the operator (|29|) after the change 
(j34j) to vanish), and the boundary condition 

lim ^ = 1 , (36) 
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i.e. so that the isotropic radius h(r) = f approaches r at spatial infinity, because 
far away from the source we assume to have euclidean geometry. This leads us to the 
definition of the function 

w(p,z) := uj [h~ x ^a/ p 2 + z 2 ^ , arctan (J- J J ? (37) 
or inversely, w such that 

u(r, 9) = w(h(r) sin 9 , h(r) cos 9) . 

Secondly, (in the spirit of Ref. 10) we define (with w(p,z)) the 5- lift of w : IR^xR — > R 
in flat K 5 , axisymmetric around the Xs-axis, by 



i 



1/2 



fv i — * uj such that £(x) = uj(xi, x 2 , x 3 , x 4 , x 5 ) := w J p = ( a^ 2 ) . : = .r- 

,*=i 

and, for every function u : M 5 — > R, the meridional cut (in direction x\) of uj 

uj i— > io such that to(p, z) := u)(p, 0, 0, 0, 2) . 

For axisymmetric functions these are isometric operations inverse to each otherP^ After 
considering the change of variable (|34|) with (|37|) in the differential operator Z ()29|) , we 
get (remember (|25|) . e A = k/j) 



where 



e A ( r )/i(r) 2 j 3 p^w + z9 2 «;| 
io»= ^ jdppw + <9^w + -<9 p w + #(r) — f^T\ J > ( 39 ) 

-A(r) A(r) 

-e 2 [-6 + 6e 3 +ri/(r)] 
lJ ' 2/i(r) ' { } 

But, through the 5-lift (J55|) . the flat Laplacian in 5 dimensions of the "lifted" function uj 
gives exactly 

5 3 
= c^u; = <9 pp w + <9 2Z w + - <9 p w , (41) 

i=i P 
first three terms in the bracket of (}3*9"j) . Furthermore, as outlined in Ref. 10, n-lift and 
meridional cut (of axisymmetric functions) leave the regularity properties and the norm 
invariant; and axisymmetric operations, like multiplication, d f (f = h(r) = (p 2 + z 2 ) 1 ^ 2 = 
($2i=i ^i 2 ) 1 ^ 2 )) an d scalar product, commute with n-lift and meridional cut. In particular, 
in the fourth term in the bracket of ()39j) the factor 

5 

pd p w + zd z w = dfUJ = '^^Xidiij . (42) 

i=i 

Therefore, substituting (jUJ) and into (03), Eq. (J2H1) in the form Z u; = -^ 2 (fi - w) 
(with Q defined from Q as it was uj from and \l/ 2 = e A 167r[eo + Po) ) now writes 

= eHr ^} r)2 {^ + g(r) E ^^ } = e A « {-167T[e (r) + po(r)][n - 0]} , 
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or, equivalently,^ 



,5 „ Ji 



ACo + H(r) E ^ X f lUJ = -lQn-f-[s (r)+ Po (r)][n - Co] . (43) 
h(r) h{r) z 

V. PROPERTIES 

With the assumptions made in Subsec. II. A for this slowly rotating configuration (starting 
from a non-rotating one as described in Subsec. IV.A.l; particularly, satisfying the energy 
condition eo + po > ), and considering only solutions of Eq. (j2H| satisfying the boundary 
and matching conditions (}3*T]) . (|3*2^) . and (J3l?j) . the following results hold 

Property (a) (Positiveness of the dragging rate) 

If the distribution of angular velocity of the fluid is non-negative (and non-trivial), then 
the dragging rate (to first order in the fluid angular velocity) is positive everywhere, 

Proof. We have seen in the former section that Eq. ()28)1 for uo is equivalent to Eq. (|4*3*|) 
for Co using the coordinate change ()34j) and the 5-lift ()38)1 . The isotropic radius f = h(r) 
is Gaussian coordinate with respect to the star's surface, f = h(R), and, thus, Co is 
at least class C 1 across this surface; therefore, from conditions on the stellar model, 
Co E W^CR 5 ) n C^M 5 ), and, hence, Eq. (02)), i.e. 

5 2 2 

LCo := ACo + J2 H ( r )j^ d ^- 16n l £ o + Po]j^y^ = - 167r N+Po]^pfi, (44) 



8=1 



is satisfied in M 5 in a generalized sense (cf. Appendix B) . Equation may be obviously 
written in divergence form 

LCo = di[a ij '(x)djU + o*(x) Co] + 6*(x) + c(x) = </(x) 

(where repeated indices denote summation over the index), with the coefficients 

a y (x) = Sij (— 1 if i — j, and = otherwise) , 
a 4 (x) = 0, 

6 l (x) = H(r)^ (Vi,je{l,...,5}), and (45) 

2 

c(x) = -167r[e (r)+po(r)]-^ (< 0) , 
and ^(x) = c(x)fi(x). (46) 

We shall consider the domain G defined by a ball in M 5 centered at the origin x = (f = 
0) and of arbitrary large radius a, 

G := B a (0) C M 5 . (47) 

Notice, whenever Cl > 0, we have, by (jjBj) , g < (because c < 0), and, hence, La) < 0, 
specifically Co e iy 1,2 (G) n C 1 (G) is a generalized supersolution relative to the operator 
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L, in (|44|) . and the domain G, (|47p. We look at the requirements for a minimum principle 
to be applied (Appendix B). The Laplacian operator is obviously strictly elliptic, and 
the coefficients (145)1 are measurable and bounded functions on G, this shows up in the 
following: the mapping r > is bounded from above and below everywhere in [0, oo[ 
(Appendix A); e + Po is also bounded, since po is bounded and po t— > £o(po) is bounded 
in any closed interval (Subsec. IV. A. 1); consequently, the coefficient c is bounded (from 
above and below); the coefficients of the first derivatives, b\ are also bounded (from 
above and below), because the function H is bounded everywhere (Appendix A) and since 
(Vi = 1, . . . , 5) Xj 2 < Y^j=i x j 2 = h(r) 2 , we have (xi/h(r)) 2 < 1. Thus, all conditions of 
a minimum principle for generalized supersolutions relative to the differential operator L 
and the domain G hold, and, as a result of the weak minimum principle (Theorem 1 in 
Appendix B), we have 

inf Co > inf Co~ , (Co~ = min(u), 0) ) . (48) 

G QG 

But, since the radius of the ball G is arbitrary, we can make it sufficiently large (a — ► oo) so 
that, by asymptotic flatness ( lim^r)-^ Co = 0, from condition (|3Tj) and lim^oo = 1), 
Co is arbitrary small at dG, following, from (|4"H|). Co > 0. Actually, the positivity is strict, 
because if u)(xo) = for some xo G G (interior point), then £>(xq) = min Co (since Co > 0), 
and, by the strong minimum principle (Theorem 2 in Appendix B), Co would be constant 
in G; in this case, Co = const. = in G (i.e. everywhere); but Co = yields, by Eq. (jUJ), 
Q = 0, or, equivalently, Q = 0, and we are assuming that Q is non-trivial. We conclude 
then io > everywhere. □ 

Property (b) 

Suppose we perturb the non-rotating configuration (in particular, with a given equation of 
state) with two (small) different distributions of angular velocity Qi and Q 2 , and integrate 
Eq. PH|) to obtain their respective solutions for the dragging rate, iO\ and u 2 , then 



^1 > ^2; ^1 ^ ^2 => 0O\> U02 . 

Proof. This follows form the linearity of Eq. ()28|) and Property (a). □ 

We are already in position to get a result about the positiveness of the difference 
Q — lo, and, hence, of the angular momentum density (122)1 . However, in order to first do 
this more specific and concrete, we shall make use of a property for the particular case of 
rigid rotation (RR), which can be found in Ref. 1, Sec. IV. 

Property RR 

For the slowly rotating configuration, 

u>(r, 0) = uj(r) 

Q{r, 9) = const. = Q > =>• < u{r) <tt in [0, R] 

{in [0, R] X [0, tt] ) w>0m[0,oo[, w'<0 in ]0, 00 [ , w'(0) = . 
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Property (c) (Positiveness of the angular momentum density) 

For the slowly rotating configuration, with a given equation of state, its dragging rate u> 
will satisfy 

co{r,e) < Q{r,6) 
if Q = Q(r, 9) (> 0) is bounded in the form 

= n < tt(r,9) < H, 

(in [0,R] x [0, 7r]) where Q is an arbitrary positive constant < Q (<^C VL crit ) (Sec. IV), 
and Q = uj(0), W solution of Eq. (j2Hl) with Q(r,9) = const. = Q, and with the same 
0-order coefficients (same starting non-rotating configuration) as the ones considered for 
our slowly rotating configuration, in particular with the same equation of state. 

Or, more generally, if (with that notation) 

u(r) < n(r,6) < H. 

Notice, u < Q means that the angular momentum density to first order in the fluid angular 
velocity, f|22|) . of this configuration (with the energy condition (fT9|) ) is > 0, vanishing on 
the axis. 

(Remarkably, the upper bound required on Q is not restrictive, because for Q continuous, 
Q is essentially bounded (g L°°) there, and fi/HfiHoo < 1.) 

Proof. We give a practical method of construction in two steps: 

1st step: Consider Q(r, 9) := Q = const. > 0, and solve the corresponding Eq. ()28|) for 
ZJ. Then, by Property RR, the solution satisfies 



u(r, 9) = ZJ(r) , (49) 
< u(r) < H in [0, R] , (50) 
and ZJ > Oin [0,oo[, uf < in ]0, oo[ , o/(0) = . (51) 

2nd step: Consider a slowly rotating configuration starting from the same non-rotating 
configuration (as in the first step) with a fluid angular velocity distribution fl(r, 9) 
such that 

U(0) =: tt < tt(r,9) < H. (52) 

Observe, we are always allowed to do this because of (joTIj) . Or, more generally, such 
that ZJ(r) < Q(r,9) < Q; notice, from ()51|) . uj is a (positive) decreasing function; in 
particular, EJ(0) > a;(r) > 0, Vr e [0, oo[. 

Now, form the second inequality in (JS*2*j) . i.e. from Q(r, 9) < Q, and, since we have the 
same starting unperturbed configuration (same 0-oder coefficients) for these both 
slowly rotating configurations, it follows, by Property (b), that their corresponding 
solutions (of Eq. (|2"Kjl ) satisfy 

u{r,9) <u7(r), (53) 
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where we have used f!49|) . On the other hand, the first inequality in (j52j) . and (foTj) 
yield 

57(r) < 57(0) := < fi(r, 0) , (54) 
and consequently, from (J53*|) and (J54*j) . 

w(r,0) < fi(r,0). 

□ 

Remark. Notice, the same argument also assures that, given a slowly rotating configura- 
tion with Q(r, 9) such that the corresponding dragging rate 5J(r, 9) < Q(r, 0), we shall have 
the same positivity result (Property (c)) for any slowly rotating configuration, starting 
from the same unperturbed configuration (in particular, with the same equation of state), 
with an angular velocity distribution Q(r, 9) such that 

fi(r, 0) := 57(r, 9) < Q(r, 6) < H(r, 0) , 

because we obtain, form the last inequality and Property (b), u(r,9) < 57(r, 0), and, 
hence, ui(r, 9) < Q(r,9). 

Series expansion. M rot . 

Before we prove next property, we first stress that, since Q and uo transform like vectors 
under rotation, Eq. (|27j) may be separated by expanding them as 



Q(r,9) = Q(r,x) ~ j^ft,(r) and (55) 

i=i 

oo 

^2^i(r)yi(x) , (56) 



u{r, U) = uj{r,x) ~ 



l=i 



with the change of variable 9 h- > x := cos0, and 

yz(x) := — — Vx G [—1,1] (9 G [0, 7r]), P; = Legendre polynomial of degree I . (57) 
ax 

Then the equation for ui takes the form 

^-[r 4 j(r) w{] + [4r 3 j» - r 2 fc(r) A,] u t = 4r 3 j"(r) ft,(r) , (58) 

with A; := 1(1 + 1) -2, I G N, I ± 0, and j and fc defined in (I25jl . 

From conditions (j3~Tj) and (}3*2j) on cu, we have the respective boundary conditions on cjj 

lim wj(r) = , (59) 

r— >oo 

^ C 1 -regular at the origin, (60) 
and, from ()33|) . the matching condition 

ui class C 1 across r = R. (61) 
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In Subsec. IV. D an explicit expression for the expansion of the rotational mass-energy 
M rot in powers or the angular velocity parameter was obtained or, using Eq. (J2H|) . 

M Iot = -- I drr 3 f [ dd sin 3 # Q(Q - uo) + 0(/i 4 ) . (62) 

Now, using the series expansions of Q and u (I55j) and and the fact that the system 
{yi}^i is orthogonal in the Hilbert space L 2 p ([— 1, 1]), with respect to the weight function 
p(x) := 1 — x 2 , x G [—1, 1], and have norm ||?/z|| p 2 = 21(1 + 1)/(2Z + 1), the integral over 6 
in (p)2"j) may be expressed as the sum 

f d0 sin 3 fl fi(r, 0) [fi(r, 5) - u(r, 9)] = V ^±il ^( r ) [ft,(r) - w,(r)] , (63) 
Jo J=1 ^ + 1 

and, consequently, the rotational mass-energy ()62|) can be expressed as a sum of integrals 
(over r) 

M - = E^TT)^ + ^ 4 )' ( 64 ) 

/■« 

with Mi := / / 2 (r) fi,(r) [n,(r) - cj,(r)] rfr , / 2 (r) := -r 3 j"(r) (> 0) . (65) 



Property (d) (Positivity and upper bound on the rotational energy M rot ) 
We consider Eq. (|58|). which can be written 

|-(r 4 j W |) -r 2 kXiUi = -4 / 2 (ft, - u;,) . 



(66) 



The main observation is that, multiplying both sides of Eq. (|66|) by uj\ and integrating from 
r = to r = oo, and taking into account that f 2 = — r 3 f = An r 4 (e +Po) k = Vr > i?, 



-^(r 4 jc^') ui-r^kXiUJj 2 



R 



dr = -4 / / 2 (ft, - w,) dr 



(note, the integral on the left hand side converges, because an asymptotically flat (JH^j) 
solution of Eq. must behave as r — > oo a;, = 0(r~ l ~ 2 ), and = 0(r _/ ~ 3 ), I > 1); 
and, after integrating once by parts the first term on the left hand side, 



4 ■ / 1°° 
r JUlVl\o 



[r 4 j ( W ;) 2 + r 2 k\i 2 ] dr = -A f u t (ft, - w,) dr . 



Now, the first term vanishes because falls off rapidly enough at r — > oo, and the second 
term (minus the integral on the left hand side) is non-positive (since j and k are always 
positive), therefore 



R 



/ 2 u;,(ft,-u;,)dr >0. 



(67) 
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We introduce now the sesquilinear (indeed bilinear) form 



R 



(u,v) f := f(r)u(r)v(r)dr u, v <E C°([0, R}), 



and the induced semi- norm \\.\\f : = ((. ,.)/) 1/2 . With this definition, M t flU} may be 
written 

Mi — (Qj, — uji)f. (68) 
We have (JHTj) . which now reads (ui, Qi — u>i)f > 0, i.e. 

(uM)/ > IMI/ , (69) 

in particular 

(flj, £*;,)/ = Kfli)/>0. (70) 
Using the Cauchy-Schwarz inequality, 

(ui,Qi)f < \\uJi\\ f \\Qi\\ f , 

which, together with (jEHJ), yields ||^||/ < ||^||/||0,||/ and, hence, \\£li\\f > \\wi\\f for 

\\u)i\\f 7^ 0, but this inequality is still true for ||co>z||/ = 0, because then, by Eq. (JH%|) . 
|| Cli\\f = 0; therefore 

IMI/ > IMI/ (71) 

holds in general. Multiplying (J7Tj) by ||Oj||/, we find, again using the Cauchy-Schwarz 
inequality, 

11^/11/ > (flljWj)/, 

but 

M l = (n h n l -ui)f = \\toi\\f 2 -(toi,wi)f, (72) 

thus showing 

^ > . (73) 
Furthermore, since (by ()70jl) (Qi,ui)f > 0, from (J72J) we also have the upper bound 

< 110,11/ , i.e. 

r R 

Mi< / / 2 (r) fi, 2 (r)dr. (74) 
/o 

The bounds (J75J) and (|71jl on M; yield respective bounds on M rot 

o < M rot < f; J|±il jf * f (r) ^ (r) dr + o(^) , 

or, writing the sum as integral over 9 (as in (|63j)). 

< M rot < i f drf(r) [ dd sin 3 fl [Q(r,9)] 2 + 0(/i 4 ) , (75) 
4 Jo Jo 

where / 2 := -r 3 / = 47r r 4 (e Q + p )e (A_I/)/2 . □ 
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VI. CONCLUDING REMARKS 



Summing up, it has been seen that relativistic stars rotating slowly and differentially, 
with a non-negative (and non-trivial) angular velocity distribution, £](x 2 ,#3) > 0), 
and satisfying the energy condition e + p > 0, have positive rate of rotational dragging 
uj > (Property (a) in Sec. V); and certain upper and lower bounds on Q assure also 
the positivity of the difference Q — uj and, hence, of the angular momentum density, 
this later vanishing on the axis, (Property (c)). We also observe that, the rotational 
mass-energy, (from Property (d)) non-negative and (as expected) "increased" by a (slow) 
angular velocity of the fluid, Q, is "decreased" by the dragging effect (over what it would 
be is this effect were neglected), i.e. is decreasing with respect to dragging rate, uj, despite 
of (as shown in Property (b)) uj being an "increasing function" of Q. 

We have given here a much simpler proof of Property (d) than in Ref. 2; this al- 
ternative proof can be even generalized outside the slow rotation limit Property (b) 
and, hence, also Property (c), however are based on the linearity of the time-angle field 
equation component to first order in the fluid angular velocity. Moreover, in the general 
differentially rotating case the rotation profile Q cannot be freely chosen, but is restricted 
by the integrability condition of the Euler equation, i.e. by Eq. (jSJ. This makes unlikely 
a generalization of Property (c) outside the slow rotation limit, other than in the form 
given in Ref. 12, Subsec. IV. B. 
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APPENDIX A: Boundedness of some functions in [0, oo[3 r 
The ratio radius - isotropic radius x(r) := 



We have 

h'{r) _ e A M/ 2 
h(r) r 

with 



(Al) 



where 



e-*> = 1 - ^ , (A 2) 
r 

m(r ).-S ^Io^(s)fds : re[0,R] 

m lr)-\ M = 4n J Q R e (s)s 2 ds : r e]R,oc[ ' [A 6) 

if we denote the stellar radius of the static model by R > 0. As we start with a (physically) 
regular (i.e. non-collapsed) static solution, we assume that 2m(r) < r (for all r e]0, R[), 
and 2M < R. 
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Integrating Eq. (jA 1|) and using Eq. (jA 2|) we get 

h(r) = h(R) exp I f , ds ) . (A 4) 

\Jr y/s(s-2m(s)) J 

(Note, the constant h(R) is determined by the asymptotic condition 

h(r) 
lim = 1 ; 

r~>oo r 

see below.) Let 

ff (r) := T -==£== Vr>0. (A 5) 



With this definition the solution, (|A 4|) . now writes 

h(r) = /i(iT) exp(#(r)) . (A 6) 

Due to the assumptions made for m, g is obviously a continuous function in the open 
interval ]0, <x>[; consequently, by (|A 6|) . h is also a continuous function there, and, in 
particular, h(r) cannot be zero in ]0, oof (unless h(R) = 0, however this would contradict 
asymptotic flatness); therefore x : r >— > x{ r ) := jjj^ is continuous in ]0, oof as well. Choose 
an e e]0, R[ and an e' &)R, oof, then x is bounded below and above in the interval [e, e'] 
(where the upper and lower bound depend on the selected e and e', of course). Let us 
now consider the intervals [0, e] and fe', oof separately: 

On fe', oof: We have 

< m(r) = M . 

Then 

- < = = . Vr G e', oo . 

s y/ s (s - 2m(s)) y/s(s - 2M) 

As in this interval r > R, we find, with Eq. (jA 5j) . 



Ry Jr s Jr y/s{s - 2M) ' \y/R+y/R-2M / 

Inserting it into Eq. (|A 6j) . yields (since exp is a monotonically increasing function) 



* r < h(r) = h(R) -7= -^=^= < — ^r. (A 7) 



R \VR + VR- 2M 



(v 7 ^ VR-2M 



Note especially that linv^oo = 4/t ^ , but, by asymptotic flatness, 

I v -R+V R—2M J 

lim^oo^ = 1; therefore h(R) = \ (V~R + VR- utj 2 > 0. Thus, h(R) > 0, 
and, from Eq. (I A 711 . 



(Vr+VR-2m) 2 r 
0< ~ WR)- X{r) -KR) <0 ° VrG[e '°° [ - (A 
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On [0, e]: We have 



< m[r) — An £o{s)s as < —Eor" : — r 
Jo 3 



2 



where i := sup rg [ fl i e (r) > 0. Now, choose e > 0, such that 1 — c r 2 > on [0, e] 
(e.g., e := (2 V ^)" 1 ). Then 

11 1 . . 

- < = < . Vr G 0, e , 

s ~ v / s {s-2m{s)) ~ sVl-cos 2 J 

and, since in this interval r < R, we find, with Eq. flA 5|) . 



In (*) = f - < -*W < f / g - m f ^ + In f 1±^^£) . 
Again, inserting it into Eq. IjA 6jl . yields 

W ^ w KR) 1 + VI - c R? ^ h(R) (1 + VI - c R 2 ) 

r > n(r) > r > r , 

i? ~ V ; ~ R 1 + VI - c r 2 - 2i? 

and, hence, 

0< ^)^ (r> ^ W (iwi-^ ) < °° Vre[M ' (A9) 

We can therefore conclude that, since IRq" = [0, e] U [e, e'] U [e', oof and x is bounded (from 
above and below) in each of these subintervals, x is bounded (from above and below) in 
K+. □ 

The function H 



We have 



and 



H(r) :-- 



— A(r) A(r) 

-e 2 [ — 6 H- 6 e 2 



2h(r) 



-A(r)/2 



rv'ir) 



1 - 



2m(r) 



2m(r) + 87rr 3 p (^) + 87rr 2 p (^) 



r — 2m(r) 



2m(r) 
r 



Thus, 



#(r) 



2/i(r) 



2m(r) 



2m(r) 



+ 87rr 2 p (^) 



2m(r) 



(A 10) 
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Now, using the Cauchy-Schwarz inequality in (jA 10|) . and the following estimates in r G 
[0, e), for some e g]0, i?[ small, (see former Subsec. in Appendix A) 



< m(r) < —r 3 



x 



< Vl-x < 1 - | Vx G [0, 1] 

< Cif < /i(r) < C2r 

< po(r) < p := sup po(r) Vr > , 

rs[0,-R] 

where the constants q (j=o,...,2) are all strictly positive (and finite), we get 

2m(r 



< 



2h(r) 
1 

2cir 

o 



6 



1 - 



2m(r) 



+ 87rr po(^) 



2m(r) 



r 



c±r, 



(All) 

(A 12) 

(A 13) 
(A 14) 



(A 15) 



with < c 3 , c 4 < oo. Therefore H is bounded in [0, e]. (Especially, due to Eq. (jA 15|) . 
H(0) = 0.) And, since, by Eq. (|A 10j) . if is also continuous in the open interval ]0, oo[ 
and lim^oo H{r) = (because lim^oo = 1), H is bounded everywhere in [0, oof. □ 



APPENDIX B: The minimum principle for generalized supersolutions 

Consider in a domain (open and connected set) G C M™ (n > 2) the differential operator 
with principal part of divergence form, defined by 

Lu = di[a.ij(x)djU + a^x) u) + bi(x) d(a + c(x) u , 

with a,ij = aji. Notice, an operator L of the general form Lu = aij(x)dijU + bi(x)diU + c(x)u 
may be written in divergence form provided its principal coefficients 5y are differentiable. 
If furthermore the constants, then even with coinciding coefficients (a^ 



c = c) and a,- = 0. Let us assume that 



1. L is strictly elliptic in G, i.e. 3 a constant A > such that A < the minimum 
eigenvalue of the principal coefficient matrix [ay(x)], 



A \y\ 2 < aij(x) yiyj Wye 



Vi 6 G; 



(B 1) 



2. ajj, a.j, 6j, and c are measurable and bounded functions in G, 



\dij\ < oo, |<2j| < oo, < oo, |c| < oo in G (i, j G {1, . . . , n}) . (B 2) 



By definition, for a function u which is only assumed to be weakly differentiable and 
such that the functions a^dju + a(a and bidiU + cu, i — 1, . . . , n are locally integrable 
(in particular, for u belonging to the Sobolev space W /1,2 (G)), u is said to satisfy Lu = g 
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in G in a generalized (or weak) sense (g also a locally integrable function in G) if it 
satisfies 

C(u,(p;G) := / {(ciijdjU + diu)di<p — (bidiU + cu)<p}dx 
Jg 

= - [ gifdx, Vy?>0 ip G C\{G) 
Jg 

(where C^(G) is the set of functions in C 1 (G) with compact support in G). 

Notice, u is generalized supersolution relative to a differential operator L and the 
domain G (i.e. satisfies Lu < in G in a generalized sense) if it satisfies £(u,<p;G) > 
0, Vy> > 9? G Cl(G). 

Theorem 1: (weak minimum principle) 

Let u G W 1,2 (G), G a bounded domain, satisfy Lu < in G in a generalized sense with 

I {op - ctdup) dx<0, \fip>0 ipe C l c {G) . (B 3) 

Jg 

and conditions ()B ljl and ()B 2j) above, 
then 

nfinu>minM~, (u~ = min(u, 0) ) . 

G 9G 



(A proof of this theorem can be found in Ref. 13, Theorem 8.1.) 
Theorem 2: (strong minimum principle) 

Let u G W 1,2 (G) PI C°(G) satisfy Lu < in G in a generalized sense, with the operator L 
satisfying conditions (|B 1)1 . fjB 2jl . and ()B 3jl . 

i/ien k cannot achieve a non-positive minimum in the interior of G, unless u = const. 

(A proof of this theorem can be found in Ref. 13, Theorem 8.19.) Note that the weak 
minimum principle, Theorem 1, for C°(G) supersolutions is a direct consequence. 



References 

1 J.B. Hartle, Ap. J. 150, 1005 (1967). 

2 J.B. Hartle, Ap. J. 161, 111 (1970). 

3 R.M. Wald, General Relativity (University of Chicago Press, Chicago, IL, 1984), pp. 
219, 220. 

4 B. Carter, Commun. Math. Phys. 17, 233 (1970). 

5 W. Kundt and M. Trimmer, Z. Phys. 192, 419 (1966). 

6 J.B. Hartle and D.H. Sharp, Ap. J. 147, 317 (1967). 

7 A. Komar, Phys. Rev. 113, 934 (1959). 



22 



8 J.M. Bardeen, in Black Holes, Les Houches 1972, edited by C. DeWitt and B.S. DeWitt 
(Gordon and Breach, New York, 1973), pp. 241-289. 

9 K.S. Thorne, in General Relativity and Cosmology, edited by R.K. Sachs (Academic 
Press, New York, 1971), pp. 237-283. 

U.M. Schaudt, Commun. Math. Phys. 190, 509 (1998). 

1 Indeed, the first order terms (in the fluid angular velocity) of the general metric, 
as given in Ref. 12, yield (omitting here the ' ~ ' symbol for all 5-lifted functions) 
H{r) = (3 dfB — 4dfU)\ h ^ , where f = h{r), given by Eq. (j3*B]) . (This can be seen 
by a straightforward calculation, using e 2U ^ r ^ = e v{ - r \ h{r) 2 e 2 ^ B{ - h ^~ u{ - h{ - r ^ = r 2 , 
and their differentiations with respect to r.) And since (from h(r) = (%2i=i x f) 1 ^ 2 ) 
dih(r) = %ih{r)~ l , we have H(r) Xih^r)' 1 = (3diB — 4diU)\ , n s , and, hence, 

Vcu) as second term in Eq. ()43j) . Also, since to first order (spheri- 
cal; K = B) e 2K N~ 1 = e 2(B ~ u ^ =r 2 h~ 2 , the coefficient of the right hand side of Eq. (gSJ) 
is actually — \l/(r) 2 e~ x ^r 2 h(r)~ 2 = — ip 2 \ ^ , in the notation of Ref. 12. 

2 M.J. Pareja, preprint gr-qc/0309073 
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